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Summary 

Two  methods  are  described  to  account  for  varying  temperature  during  creep. 
Both  employ  the  modified  superposition  principle.  One  uses  a  reduced  time  involv¬ 
ing  a  shift  factor  which  is  a  function  of  both  stress  and  temperature  history. 

The  other  considers  the  strain  to  be  a  function  of  the  current  values  of  stress 
and  temperature. 

Experiments  on  polyurethane  include  constant  stress  creep  and  recovery  at 
several  temperatures  in  the  nonlinear  range  and  an  experiment  in  which  the  stress 
was  held  constant  while  the  temperature  increased  at  a  constant  rate,  then  the 
stress  was  removed  and  the  temperature  decreased  at  a  constant  rate.  The  strain 
in  this  experiment  was  predicted  by  means  of  the  theories  from  the  results  of  the 
constant  temperature  creep  tests. 

The  strains  in  the  constant  temperature  creep  and  recovery  tests  were 
described  by  means  of  the  multiple  integral  representation  and  the  modified  super¬ 
position  principle.  Most  of  the  nonlinearity  and  temperature  effect  were  found 
in  the  coefficient  of  the  time-dependent  term. 
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Introduction 


In  a  recent  paper  [1]  the  authors  described  the  results  of  creep  experi¬ 
ments  in  the  nonlinear  range  on  polyurethane  in  torsion  at  various  temperatures. 

The  results  showed  that  the  time  dependence  could  be  described  adequately  by  a 
power  function  whose  exponent  was  independent  of  stress  and  temperature.  In  the 
present  paper  tension  creep  experiments  in  the  nonlinear  range  performed  on  the 
same  specimen  of  polyurethane  as  part  of  the  same  series  of  experiments  are  re¬ 
ported  for  several  temperatures .  Constitutive  equations  describing  these  results 
are  employed  to  predict  the  behavior  under  a  constant  rate  of  temperature  rise  (and 
fall)  during  creep  (and  recovery)  at  constant  stress. 

Previous  investigations  of  nonlinear  viscoelasticity  have  mostly  been  con¬ 
fined  to  room  temperature.  However,  Bernstein,  Kearsly  and  Zapas  [2]  and  Lianis 
[3]  have  extended  the  concept  of  a  thermorheologically  simple  viscoelastic  mate¬ 
rial  to  the  nonlinear  range.  Earlier  investigations  by  the  authors  [4  to  9]  on 
poly  (vinyl  chloride)  and  polyurethane  at  room  temperature  under  combined  tension 
and  torsion  showed  that  the  creep  strain  could  be  separated  into  time-dependent 
and  time-independent  parts .  The  latter  could  be  expressed  as  a  power  function 
whose  exponent  was  independent  of  stress  state.  The  time-independent  term  was 
nearly  linear,  but  the  time -dependent  coefficient  was  strongly  nonlinear.  It  was 
found  that  the  nonlinear  behavior  could  be  adequately  described  by  the  first  three 
terms  of  a  multiple  integral  series.  For  varying  loads  either  the  modified  super¬ 
position  method  [5,10]  for  the  mixed  time  parameters  in  the  kernel  functions  or 
the  product  form  [6,10]  gave  satisfactory  results. 


Numbers  in  brackets  identify  references  listed  at  the  end  of  this  paper. 


2 


Material  and  Specimen 

The  material  used  in  the  experiments  was  the  same  as  employed  in  previous 
work  [8,9,11,12],  It  was  a  full  density  solid  polyurethane  having  a  specific 
gravity  of  1.25.  This  material,  identified  as  XR6-77,  was  produced  by  CPR  Divi¬ 
sion,  the  Upjohn  Company,  and  Supplied  by  the  Lawrence  Radiation  Laboratory.  The 
material  was  prepared  from  the  prepolymer  and  combined  with  toluene  diisocyanate. 
The  polyester  resin  was  derived  from  a  dicarboxylic  acid  and  a  triol.  The  carbon 
content  was  59  per  cent,  hydrogen  6.5  per  cent,  nitrogen  7.5  per  cent,  chlorine 
less  than  0.3  per  cent,  ash  0.05  per  cent,  and  the  balance  oxygen.  The  material 
was  cured  at  a  temperature  of  250°F  for  two  hours. 

The  specimen  was  a  tubular  type  with  enlarged  threaded  ends  machined  from 
a  solid  rod.  The  average  outside  diameter  was  0.9971  in.,  with  a  maximum  devia¬ 
tion  of  0.0003  in.,  the  average  wall  thickness  was  0.05928  in.,  with  a  maximum 
deviation  of  0.0005  in.,  and  the  gage  length  was  4.00  in.  All  experiments  were 
performed  on  the  same  specimen. 

Experimental  Apparatus 

The  testing  machine  was  designed  to  test  tubular  specimens  under  combined 
tension,  torsion,  and  internal  pressure  [13].  It  has  been  used  to  investigate 
creep  and  stress  relaxation  behavior  under  various  loading  conditions  for  plastics 
[4-12]  and  metals  [14-16]  and  plasticity  [17]. 

The  loading  and  measuring  devices  are  as  follows.  Tension  and  torsion 
were  produced  by  dead  weights  or  servo  control  acting  through  levers  and  pulleys. 
Dead  weights  were  used  in  the  present  work.  The  tensile  strain  measurement  was 
accomplished  by  using  pairs  of  rods  attached  to  the  upper  and  lower  gage  points 
of  the  test  specimen  to  transmit  the  motion  of  the  gage  points  below  the  furnace. 
The  lower  ends  of  one  pair  of  the  rods  were  attached  to  the  core  and  the  ends 
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of  the  other  pair  were  attached  to  the  coil  of  a  differential  transformer.  Rela¬ 
tive  motion  of  the  upper  and  lower  gage  points  during  the  test  caused  an  output 

of  the  differential  transformer  which  was  directly  related  to  the  strain.  The 

—6 

sensitivity  of  the  strain  measurement  was  2  x  10  in. /in. 

In  order  to  perform  creep  tests  of  polyurethane  above  room  temperature, 
the  following  heating  devices  were  added.  The  specimen  was  heated  in  an  electric 
split-tube  furnace  which  was  maintained  at  constant  temperature  by  a  C.  N.  S. 
Instruments'  Sirect  Mark  III  proportional  temperature  controller.  Due  to  the  poor 
conductivity  of  polyurethane  a  second  heating  source  inside  the  specimen  tube  was 
used  to  achieve  a  more  uniform  temperature  distribution.  Two  small  resistor  heat¬ 
ing  cores  were  attached  at  both  ends  of  a  copper  bar  and  inserted  inside  the 
specimen.  It  was  found  that  the  smaller  the  spacing  between  the  copper  bar  and 
the  inner  surface  of  the  specimen  the  better  the  temperature  was  controlled.  One 
thermocouple  was  embedded  in  each  end  of  the  copper  bar.  These  served  as  the  heat 
sensing  device  in  the  control  system.  Temperature  in  the  copper  bar  was  controlled 
by  two  Research  Incorporated  Thermae  Model  TC5192  temperature  controllers  operating 
on  the  two  heaters . 

During  each  of  the  tests,  temperature  was  measured  at  six  positions  by 
means  of  chromel-alumel  thermocouples  with  fused  hot  junctions.  Two  of  these 
thermocouples  were  cemented  on  the  inner  surface  of  the  specimen  at  positions 
slightly  beyond  the  gage  length.  Four  thermocouples  were  attached  with  spring 
clips  on  the  outer  surface  at  equal  distances  within  the  gage  length  of  the  speci¬ 
men.  It  was  observed  that  the  temperature  could  be  maintained  constant  throughout 
the  testing  period  within  ±  1/2°F.  The  temperature  variation  along  the  gage  length 
was  uniform  at  75°F  and  102°F,  but  increased  to  2°F  on  the  outside  and  7°F  on  the 
on  the  inside  of  the  specimen  at  150°F. 
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Experimental  Procedure  and  Results 

The  rate  of  heating  the  specimen  was  controlled  to  about  2°F  per  minute. 
After  reaching  the  test  temperature  the  specimen  was  kept  at  this  temperature  for 
one  hour  before  the  load  was  applied.  Pilot  tests  indicated  that  a  soaking  period 
from  0.5  hr.  to  5.0  hr.  at  the  test  temperature  (75°F  to  160°F)  did  not  show  a 
significant  difference  in  the  creep  response  under  the  same  loading  conditions. 
Except  for  one  experiment  the  temperature  was  maintained  at  a  constant  value 
throughout  the  testing  period.  Some  adjustments  of  the  set  points  of  the  con¬ 
trollers  were  required,  however. 

In  performing  the  constant  stress  creep  tests,  the  specimen  was  loaded  for 
a  one-hour  period  for  each  test.  After  unloading  at  the  end  of  one  hour,  the 
specimen  was  left  to  recover  at  zero  stress  at  the  same  test  temperature  for  about 
one  hour.  After  that,  heat  was  gradually  turned  down  until  it  reached  room  tem¬ 
perature,  and  the  specimen  was  left  to  recover  at  room  temperature  until  the  strain 

~6 

had  returned  to  less  than  40  x  io  in. /in.  or  appeared  to  have  stabilized  before 

proceeding  to  the  next  test.  Tests  69;  46;  48;  84;  91;  and  97  followed  tests  whose 

recovery  stabilized  at  =  102  ;  =  78  ;  =  95  ;  =  120  and  =  00 

~6 

=  -84  ;  and  =  63  x  10  in. /in.,  respectively.  The  zero  strain  for  each 

test  was  taken  to  be  the  strain  reading  at  temperature  just  prior  to  loading. 

Thermal  expansion  has  two  effects  on  tensile  strain:  it  shows  up  as  an 
apparent  creep  strain  during  temperature  changes ,  and  it  alters  the  gage  length  of 
the  strain  measuring  instrument.  The  former  was  taken  into  account  in  the  con¬ 
stant  rate  of  temperature  rise  test,  but  the  latter  was  ignored  in  all  tests.  The 

gage  length  was  fixed  at  75°F  and  no  correction  was  made  for  change  in  gage  length 

—6 

with  temperature.  The  thermal  expansion  was  found  to  be  35.2  x  10  in./in./°F. 

Thus  at  160°F  the  gage  length  increased  about  0.3  per  cent. 
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The  test  program  of  the  present  investigation  is  shown  in  Table  I .  The 
tensile  strain  versus  time  curves  from  these  tests  are  shown  in  Fig.  1  and  2.  A 
constant  stress  creep  test  was  also  conducted  in  which  the  temperature  increased 
at  a  nearly  constant  rate,  then  decreased  at  a  nearly  constant  rate  while  the 
stress  was  removed  when  the  temperature  started  to  decrease.  The  strain  versus 
time,  together  with  the  programmed  and  measured  temperature,  are  shovm  in  Fig.  3. 

In  previous  work  on  the  same  material  at  room  [8,9,11]  and  elevated  [1] 
temperature,  it  was  found  possible  to  describe  the  results  of  creep  strain  versus 
time  under  constant  stress  by  a  power  law  with  a  constant  exponent  n  , 


A  +  .n 

♦  'n* 


CD 


where  and  are  functions  of  stress  and  n  is  a  constant.  From  the 

* 

results  shown  in  Fig.  1  and  2,  it  seems  that  (1)  can  be  used  to  describe  the 


strain  versus  time  relationship  of  the  present  tests  also.  The  effect  of  tempera¬ 
ture  on  e°  ,  e+  ,  n  in  (1)  was  determined  by  first  obtaining  the  best  fit  of 


each  data  set  to  (1)  by  a  least  squares  method.  As  in  the  companion  tests  in  tor¬ 


sion  it  was  found  that  there  was  considerable  scatter  in  n  but  no  definite  trend 
as  a  function  of  either  stress  or  temperature.  showed  scatter  and  about  the 

same  modest  temperature  effect  observed  in  torsion  [1];  e°^/ a  increased  from 
about  0.19  to  0.23  from  75°F  to  160°F.  However,  showed  a  strong  effect; 

it  more  than  doubled  over  the  same  temperature  range,  the  same  as  in  torsion  and 
for  the  same  reason  [1] . 


In  view  of  the  scatter  in  n  resulting  from  the  individual  determin¬ 
ations  the  value  n  =  0.143  used  in  earlier  work  [8,11]  with  this  material 
at  room  temperature  and  in  the  companion  work  in  torsion  at  elevated  tempera¬ 
ture  [1]  was  chosen  as  a  common  value  for  all  stresses  and  temperatures. 

S 

Numbers  in  parentheses  identify  equations. 


was  taken  to  be  independent  of 
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Also  consistent  with  the  work  on  torsion  [1] 
temperature  in  view  of  the  much  smaller  apparent  effect  of  temperature  on 
than  on  .  Employing  these  restrictions  values  of  and  were  re¬ 

calculated  and  are  shown  in  Table  I.  In  the  analysis  in  following  sections  equa¬ 
tion  (1)  with  the  values  of  and  n  in  Table  I  were  employed  rather 

than  the  experimental  data. 

Discussion  of  Results 

Several  experiments  were  repeated  one  or  more  times  at  intervals  during 
the  series  of  118  creep  tests  under  tension,  torsion  and  combined  tension  and  tor¬ 
sion.  These  repeated  experiments  were  used  to  indicate  how  much  the  creep  behavior 
was  changing  as  a  result  of  the  history  of  stress  and  temperature  on  the  single 
specimen  employed.  As  described  in  [1]  the  effect  of  history  was  greater  in  tor¬ 
sion  than  tension.  The  results  of  the  repeated  tests  in  tension  are  summarized 
in  Table  II .  Examination  of  all  the  data  showed  that  there  were  changes  which 
generally  resulted  in  increased  strains.  The  greatest  changes  occurred  during  the 
early  tests  and  again  when  the  largest  stresses  were  applied.  Possible  causes  were 
discussed  in  [1].  There  was,  however,  an  extensive  period  in  which  changes  were 
small.  Compare,  for  example,  tests  39  and  76  and  tests  41  and  67  in  Table  II. 

In  the  analysis  of  the  effect  of  temperature  the  data  from  the  first  20 
experiments  were  ignored  because  of  the  changes  which  were  occurring.  Also  some 
of  the  last  tests  were  omitted  and  a  few  others  whose  strains  were  so  out  of  line 
as  to  suggest  errors  of  observation.  The  experiments  employed  in  the  analysis  were 
undoubtedly  influenced  to  some  extent  by  the  prior  history  of  stress  and  tempera¬ 
ture.  These  history  effects  were  ignored  in  the  analysis. 
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Analysis 


Nonlinear  Constitutive  Equation 

A  nonlinear  constitutive  equation  in  the  multiple  integral  form  has  been 
described  by  several  investigators;  see  for  example  [4],  For  uniaxial  stress 
=  a  the  following  equation  can  be  obtained  for  the  creep  strain  =  e  , 

ft 


e 


R(t-51;T)a(C1)d?1 


o 


+||  M(t-51,t-C2;T)d(?1)d(C2)d51dC2 
o  o 

+  j  I  j  N(t-S1,t-52,t-53;T) 0 (ei)d(52)d(C3)dC1d52dC3  ,  (2) 

o  o  o 

where  t  is  time,  T  is  temperature  £  an  arbitrary  time  prior  to  t  and 

6(0  =  3 o(£)/3£  .  In  this  equation  only  terms  from  the  first  three  orders  have 

been  retained.  The  kernel  functions  R,M,N  of  (2)  depend  on  temperature 

(under  different  isothermal  conditions)  and  time  arguments  as  indicated.  Under 

& 

constant  stress  a(t)  =  aH(t)  ,  where  H(t)  has  the  value  1  when  t  Z  0  and 
0  when  t  <  0  .  Thus  under  isothermal  conditions  and  constant  stress  (2) 
becomes 

e(t,T)  =  R(t,T)cr  +  M(t,T)o2  +  N(t,T)a3  .  (3) 


Determination  of  Kernel  Functions 

Comparison  of  (3)  with  (1),  representing  the  test  data,  suggests  that  the 
kernel  functions  R  ,  M  ,  N  may  have  the  following  forms 


H(t)  is  the  Heaviside  unit  function. 
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R(t,T)  =  ax  +  b^TH” 

M(t,T)  =  a2  +  b2(T)tn 

N(t,T)  =  a3  +  b3(T)tn  ,  (4) 


where  n  =  0.143  ,  a^  ,  a2  ,  a3  are  independent  of  temperature,  and  b^  ,  b2  , 
bg  are  functions  of  temperature  as  suggested  by  the  test  results.  Inserting 
(4)  into  (3)  and  comparing  with  (1)  results  in 


o  A  2  A  3 

e  =  a^  +  a2o  +  a3<j 


e+  =  b1(T)a  +  b2(T)a2  +  b3(T)o3 


( 5a) 

(5b) 


From  the  results  in  Table  I,  the  time- independent  and  the  time- dependent  coeffi¬ 
cients  (5a)  and  (5b),  respectively,  were  determined  as  follows. 

The  data  of  e  for  each  stress  as  given  in  Table  I  shows  that  e  can 
o  o 

be  approximated  as  a  linear  function  of  stress.  Therefore 


e  =  a  0  , 

o  l 

a^  =  0.195  per  cent/ksi.  ,  a2  =  a3  =  0  .  (6) 

The  time-dependent  part  of  the  strain  e+  from  Table  I  versus  temperature 
T  was  plotted  for  each  stress  level  as  shown  in  Fig.  4.  Then  e+/o  versus  o 
for  each  of  several  temperatures  (75°F,  100° F,  120° F  and  140° F)  was  plotted  as 
shown  in  Fig.  5.  In  Fig.  5  values  of  e+  for  temperatures  100°F,  120°F  and 
140°F  at  all  stress  levels  were  obtained  from  the  lines  of  Fig.  4. 

From  (5b) 

ef(o,T.)  , 

- =  b^)  ♦  b2(Tl)o  ♦  bjCT^o  , 

where  T^  =  T  -  75°F  .  The  reference  temperature  was  chosen  as  75°F  because 


(7a) 
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T  =  75°F  was  the  lowest  temperature  used  in  this  study.  In  Fig.  5,  curves  of 
t+/o  versus  a  at  different  temperatures  are  parallel  straight  lines.  Thus 
the  curve  for  different  temperatures  can  be  obtained  by  a  vertical  shift  of  the 
e+/a  versus  o  curve  for  =  0  .  This  implies  that  only  b^  depends  on  tem¬ 
perature,  b^  and  b^  are  independent  of  temperature.  Therefore  (7a)  can  be 
rewritten  as 

e+(cr,T  )  2 

- —  =  b.  +  b0o  +  b.o  +  F(T  )  ,  (7b) 

a  l  ?  6  l 

where  b^  ,  b^  ,  bg  are  independent  of  temperature,  F(T^)  is  the  vertical 
shift,  and  F(0)  =  0  .  b^  ,  b2  ,  b^  were  determined  graphically  from  e+/o 
versus  c  ,  in  Fig.  5,  at  T^  =  0  with  the  result 

b^  =  0.017,  per  cent/ksi 

2 

b^  =  0.004,  per  cent/(ksi) 

b3  =  0  .  (7c) 

F(T^)  was  determined  from  the  amount  of  vertical  shift  of  e+/a  versus  temper¬ 
ature  as  shown  in  Fig.  6.  Thus  F(T^)  was  approximated  by 

F(T  )  =  0  T  +  02T^  ,  (7d) 

where  0^  =  0.00038  per  cent/ksi  °F  and  02  =  0.000,0045  per  cent/ksi(°F)2  . 

From  (7a),  (7b),  (7c),  (7d)  the  time-dependent  part  of  the  strain  can  be  repre¬ 
sented  by 

e*(o,T)  =  (bx  +  0^  +  02T^)o  +  b2a2  .  (8) 

Thus  the  constitutive  relation  of  creep  strain  under  constant  stress, 
constant  temperature  as  characterized  in  this  section  may  be  summarized  as  fol¬ 
lows  by  introducing  (4)  through  (8)  into  (3) 
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or 


e  =  +  [(b  +  61T1  +  62T2)o  +  b2a2^tn  ’ 


e  =  0.195a  +  [(0.017  +  0.00038^  +  0.000,0045T  )o 


+  0.004a2]t°‘143 


(9) 


where  e  is  in  per  cent,  a  in  ksi,  and  t  in  hours.  Creep  strains  were  cal¬ 
culated  from  (9)  for  each  stress  and  temperature  and  are  shown  in  Fig.  1  and  2. 
Recovery  at  Constant  Temperature 

The  modified  superposition  method  [5]  for  multi-step  load  histories  can 
be  expressed  as  follows  for  tensile  loading  a 

P 

e  (t)  =  l  [f(o.,t-t.)  -  f(o.  ,t-t.)]  ,  t  >  t  (10) 

J.X  j=q  *  1“X  X  p 

where  for  the  present  experiments 

f(o,t,T)  =  R(t,T)o  +  H(t,T)o2  +  N(t,T)a3  ,  (11) 

which  becomes  (9)  when  numerical  values  are  supplied.  Thus  the  recovery  at  zero 
stress  after  a  period  t^  =  1  hr.  at  constant  stress  oq  is  from  (10)  for  iso¬ 
thermal  conditions 

e .,(t)  =  f(a  ,t)  -  f(a  ,t-t  )  ,  (12) 

11  o  o  1 

where  f  has  the  form  of  (9)  and  the  a^  term  vanishes. 

The  recovery  strains  following  constant  stress  creep  were  predicted  from 
(9)  by  means  of  (12).  The  results  for  some  of  the  tests  are  shown  in  Fig.  7. 
During  recovery  the  contribution  of  the  nonlinear  terms  is  very  small.  Agreement 
between  the  theory  and  test  data  for  tension  (Fig.  7)  is  not  as  good  as  for  tor¬ 
sion  [1].  With  the  exception  of  one  test  the  actual  strains  during  recovery  were 
smaller  than  the  strains  predicted  from  the  tension  tests. 
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Creep  under  Continuously  Varying  Temperature 

Method  A 

In  the  analysis  of  the  creep  behavior  of  linear  viscoelastic  material  at 
different  but  constant  temperatures,  the  effects  of  temperature  and  time  for  cer¬ 
tain  classes  of  material  can  be  combined  into  a  single  parameter  through  the 
"time -temperature  superposition  principle"  which  may  be  stated  as  follows: 


j(t,o  =  j(t  ,t) 

o 

(13) 

C  =  t/aT(T) 

(14) 

where  J  is  the  creep  compliance,  5  is  the  reduced  time,  a^,  is  the  tempera¬ 
ture  shift  factor,  Tq  is  the  reference  temperature  and  t  is  the  real  time. 
Thus,  the  determination  of  a^,  as  a  function  of  temperature  T  will  provide  the 
necessary  information  for  determination  of  the  reduced  time.  For  tests  involving 
varying  temperature  during  loading,  the  definition  given  by  Morland  and  Lee  [18] 
for  reduced  time  is 


(15) 


o 

where  t  is  the  current  time,  C  is  any  prior  time,  and  t  =  0  corresponds  to 
the  first  application  of  stress. 

Aside  from  the  work  reported  in  [2]  and  [3],  no  similar  method  is  avail¬ 
able  to  describe  nonlinear  viscoelastic  behavior  under  varying  temperature  during 
loading.  It  is  intended  in  this  section  to  explore  a  possible  way  of  describing 
the  effect  of  nonlinear  creep  behavior  under  varying  temperature  by  generalizing 
the  constitutive  relation  obtained  through  creep  tests  under  different  but  con¬ 
stant  temperatures. 

Consider  a  general  constitutive  relation  under  constant  stresses  and  tem¬ 


peratures,  as  follows: 


12 


e  =  J  a  +  J.  a 
o  t 


(16) 


where 


J  =  F  (T)  +  F0(T)o  +  T,(T)o‘ 
o  1  ^  O 


Jt  =  [f|(T)  +  F^Do  +  Fg(T)o2]G(t) 


(16a) 

(16b) 


In  (16)  it  is  assumed  that  the  total  creep  strain  can  be  separated  into  a  time- 
independent  part  Jq0  and  a  time-dependent  part  •  Equation  (9)  can  be  con¬ 

sidered  as  a  special  case  of  (16).  Rearrange  (16b)  such  that 

Jt  =  [“l  +  a2°  +  a3°2  +  B1(V  +  62(Tl)a  +  ^3(T1)°2]G(t)  •  (17) 

Rearranging  by  factoring  out  the  temperature- independent  terms  yields 


Jt  =  [a]L  +  o2a  +  a3a  ]  f(T1,a)G(t) 


(18) 


where 


f(T1>a)  5 


1  + 


31(T1)  +  B2(Tl)a  +  63(T1)0 

•  2 

(a^  +  a^a  +  a^a  ) 


and  a1  ,  o2  ,  ag  are  independent  of  temperature,  while  8^^  ,  B2  ,  83  are 
functions  of  temperature  =  T-Tq  ,  and  6^  =  8g  =  3  3  =  0  at  T  =  Tq  or  =  0 

f(T1,c)G(t)  in  (18)  can  be  combined  into  G (?)  by  defining  a  reduced 

time  C  , 


C  = 


aT(T1,o)  ’ 


(19) 


where  a^T^o)  is  a  shift  factor  depending  on  both  temperature  and  stress. 
Therefore , 

e  =  J  a  +  [a.o  +  a_a  +  a_a  ]G(^)  .  (20) 

o  1  2  3 

Equation  (20)  can  be  interpreted  as  a  time-temperature  superposition  principle 
[20]  provided  that  the  shift  factor  aT  depends  on  both  temperature  and  stress. 
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For  continuously  varying  temperature  and/or  stress,  (18)  and  (19)  can 
be  used  to  redefine  the  reduced  time  ?  as  follows: 


f*  _ ds 

J  aT[T1(s), 


o(t)] 


where  s  is  time  prior  to  the  current  time  t  and  t  =  0  at  the  time  stress 
is  first  applied. 

For  describing  creep  behavior  under  an  arbitrary  stress  history,  it  has 
been  found  that  the  concept  of  a  "modified  superposition  principle,"  proposed  by 
Findley  and  Lai  [5]  in  summation  form  and  employed  in  integral  form  by  Pipkin  and 
Rogers  [19],  yields  satisfactory  agreement  with  experiment. 

The  modified  superposition  method  considers  that  the  strain  following  a 
change  in  stress  at  time  t^  is  the  sum  of:  (a)  the  strain  that  would  have 
occurred  had  the  stress  not  been  changed,  (b)  minus  the  strain  that  would  have 
occurred  had  the  prior  stress  been  applied  at  the  time  t^  to  a  previously 
unstressed  material,  (c)  plus  the  strain  that  would  have  occurred  had  the  current 
stress  been  applied  at  time  t^  to  the  previously  unstressed  material. 


In  integral  form  the  modified  superposition  principle  can  be  written  as 


follows 


./^  _  f*  9Q[a(C),t-C,T]  d oU)  „ 

(t)  "  J  - a^T V - d T~dK  * 


Introducing  (20)  into  (22)  for  Q[a(0,t-C,T]  provides  a  means  of  account¬ 
ing  for  creep  behavior  under  continuously  varying  stress  and  temperature  as 


follows : 


:(t )  =  J  o(t)  +  f?  ^  (M^>  +  a.oU’)2  +  a-ouVlGuV)}  ds' 

°  J  3a  1  1  2  3  3r 


(23) 
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where 


T  ds 

aTCT1(s) ,o(t] 


C  =  \  a  TOOVoTO]  •  (25) 

o  1 

For  T  =  Tq  ,  a  constant, (23)  reduces  to  the  form  originally  proposed  for 
the  modified  superposition  method  as  follows: 

►  "£  .  . 

s(t)  =  JQ0(t)  +  {[0^0(5)  +  a 2o2(0  +  a3o3(£)]G(t-0}  ^  -  dg  » 


since  as  defined  in  (18)  for  T  =  T  ,  f(T,o)  =  1  ,  £  =  t  ,  and  C  =  £  . 

o 

For  T  >  Tq  but  constant,  and  for  varying  stress  a  =  a(t)  ,  the  strain 
may  be  described  by  (23)  provided  that 


ft  , 

ds 

a-rCv 


*tl  1  ’ 


. 


j 


aT[T1,o(e>:  • 


For  varying  temperature  but  constant  stress  0  =  0QH(t)  .  Thus  (23) 


reduces  to 


where 


2  3 

e(t)  =  J  0  +  [a,0  +  a _o  + a-0  ]G(s)  , 

OO  lO  zO  OO 


x  l'  y ’  o‘ 


and  c  is  not  involved. 


H(t)  is  the  Heaviside  function.  It  has  the  value  one  when  t  >,  0  and  zero 
when  t  <  0  . 


Equation  (30)  was  used  to  predict  the  creep  strain  (total  strain  minus 


strain  due  to  thermal  expansion)  for  the  test  program  shown  in  Fig.  3,  as 
follows.  In  this  program  the  stress  was  constant  o  =  4000  psi  ,  and  the  tem¬ 
perature  T^(t)  =  T  -  75°F  =  60t  where  t  is  time  in  hours. 

Equation  (9)  may  be  expressed  in  the  form  of  (18)  as  follows 


e(t)  =  a^  +  (b^o  +  b2a  ) 


1  + 


01T1  +  Vl 

bl  +  b2° 


Comparing  (31)  with  (18)  shows  that 


G(t)  =  tn 


f(Tl5o)  = 


1  + 


01T1  +  02T1 

bl  +  V 


Therefore 


GU)  =  Cn  =  f— )"  =  f(T  ,cr)G(t)  = 
aT  1 


1  + 


Vl  +  Vl 

bl  +  b2° 


Hence  (31)  becomes 


e(t)  =  a^o  +  (b^o  +  b2o2)?n 


(32 )  implies  that 


1  + 


Vl  +  e2Tl 

b^  +  b"2o  1  ’ 


or 


a^T^o] 


1  + 


61T1  +  02T1 
bl  +  b2°“ 


0 


For  varying  temperature  in  which  =  60t 
time  c  is  defined  from  (34)  in  (30)  as 


and  o=o  =  4  ksi  the  reduced 
o 
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rt 

C  =  1  + 

o  L 

Equation  (35)  was  evaluated  numerically  by  using  Simpson's  rule.  Then 
the  creep  strain  was  calculated  by  using  (35)  in  (33)  with  values  of  a^  ,  b^  , 
b2  ’  0i  ’  ®2  ’  n  &iven  in  (9).  The  total  strain  eT(t)  is 

eT(t)  =  aT^  +  e(t)  ,  (36) 

where  a  =  35.2  x  io  6  in./in./°F  is  the  observed  value  of  the  coefficient  of 
thermal  expansion  determined  by  measuring  the  thermally  induced  strain  at  zero 
load. 

In  (36)  it  was  assumed  that  the  thermal  expansion  coefficient  was  inde¬ 
pendent  of  the  stress. 

The  values  of  e(t)  and  eT(t)  predicted  from  (33,35)  were  calculated 
and  are  shown  in  Fig.  3.  The  strains  were  also  evaluated  numerically  from  the 
actual  temperature  instead  of  T^  =  60t  as  shown  in  Fig.  3.  For  this  calculation 
the  temperature  change  was  approximated  as  steps. 

Method  B 

A  more  direct  means  of  describing  the  effect  of  temperature  (both  constant 
temperature  and  varying  temperature)  on  the  time-dependent  part  of  the  strains 
(16b)  is  as  follows.  Consider  that  temperature  only  affects  the  coefficients 

tit 

F  (T)  ,  F  (T)  and  F„(T)  of  (16b).  Employing  the  modified  superposition  method 
X  2  3 

the  creep  behavior  under  varying  load  and  temperature  can  be  described  from  (16) 
and  the  modified  superposition  principle  (22)  as  follows 

c(t)  =  J  o(t)  +  P  ^-{JtC(t-C),o(0,T(t)3a(C)}  ^|^-d£  .  (37) 

o 


2lk 

91(60s)  +  62(60s)  n 
b1  +  b2(4) 


where 


JtC(t-?),a(5),T(t)]  =  [o1  +  a2a  +  a3c2]f  (^(t)  ,o)G(t-  O 
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For  the  present  material  the  temperature  and  stress  function  in  J  is 
given  by  the  function  in  square  brackets  in  (9).  In  the  test  program  shown  in 
Fig.  3  a  =  4000  psi  ,  and  T  -  75  =  T^  =  60t  .  Since  the  stress  was  constant, 
integration  of  (37)  by  means  of  the  Heaviside  function  and  making  use  of  (9) 
yielded  the  following  prediction  of  strains  during  this  part  of  the  experiment: 

e(t)  =  a^o  +  {[b1  +  0  (60t)  +  e2(60t)2]o  +  b2o2}tn  .  (38) 

The  corresponding  total  strain  is 

e^,(t)  =  a(60t)  +  e(t)  .  (39) 

The  results  calculated  from  (38,39)  using  the  values  of  a^  »  ^  *  ®i  ’  ’ 

n  given  in  (9)  are  also  plotted  on  Fig.  3. 

Also  shown  in  Fig.  3  are  the  strains  calculated  from  (38)  and  (39)  using 
the  actual  temperature  instead  of  T^  =  60t  . 

Figure  3  shows  that  (38)  and  (39)  (Method  B)  are  closer  to  the  experimental 
data  than  (33-36)  (Method  A).  Also  shown  by  square  points  in  Fig.  3  are  the 
strains  at  the  indicated  time  from  creep  tests  at  a  constant  stress  of  4000  psi 
and  the  constant  temperature  indicated.  Comparing  these  with  the  data  for  constant 
stress  but  increasing  temperature  it  is  observed  that  the  four  points  nearly  coin¬ 
cide  with  Method  B,  as  would  be  expected  from  the  nature  of  Method  B,  but  they 
lie  below  and  nearly  parallel  to  the  increasing  temperature  test  instead  of  above 
as  would  be  expected.  This  fact  may  result  from  differences  in  temperature  dis¬ 
tribution  between  the  two  types  of  tests.  The  prediction  of  Method  A  lies  still 
farther  below  the  test  data  for  the  increasing  temperature  test. 

Recovery  under  Continuously  Varying  Temperature 

After  the  rising  temperature  reached  the  programmed  value  of  135° F  at  one 
hour  the  temperature  was  caused  to  decrease  gradually.  Also,  the  load  was  re¬ 
moved  at  one  hour.  The  temperature  and  strain  during  this  period  of  recovery  at 
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decreasing  temperature  were  as  shown  in  Fig.  8.  The  strain  corresponding  to  the 
recovery  at  decreasing  temperature  was  computed  for  both  Methods  A  and  B  by  using 
the  modified  superposition  principle  to  account  for  the  abrupt  removal  of  the 
stress  at  one  hour.  That  is,  the  actual  history  was  considered  to  be  equivalent  to 
(a)  creep  at  constant  stress  through  the  entire  period  but  at  the  actual  variation 
in  temperature,  minus  (b)  the  creep  at  the  same  stress  applied  at  one  hour  with  the 
decreasing  temperature  that  actually  occurred  from  one  hour  on.  Since  the  tempera¬ 
ture  variation  did  not  occur  at  a  constant  rate,  especially  during  falling  tempera¬ 
ture,  the  temperature  variation  was  approximated  by  a  series  of  appropriate  steps 
of  constant  temperature.  This  method  accounted  for  the  small  variations  from  con¬ 
stant  rate  of  temperature  rise  shown  in  Fig.  3.  The  results  for  the  rising  portion 
by  the  step  method  were  nearly  the  same  as  for  T^  =  60t  as  shown  in  Fig.  3. 


Computation  of  the  recovery  strain  by  Method  A  using  the  modified  super¬ 
position  principle  was  as  follows: 


e(t)  =  (be  +  bo  )[G(?)  -  G (?  -  £,)]  , 

i  o  /  o  1 

t  >  1  hr. 


where  oq  =  4  ksi  ,  GiO  =  and  G(c  -  c  )  =  (c  -  ,  and 


ft 


C  = 


1  + 


6^(5)  +  e2T^(s)”11 


O  *“ 


?1  = 


o  *- 


1  + 


bl  +  Vo 


HjT^s)  +  e2q(s) 


bi  +  Vo 


n 


ds 


ds 


(40) 


(41) 


(42) 


where  t  =  0  was  taken  to  be  the  time  at  which  the  first  stress  was  applied,  and 
(41,42)  were  evaluated  by  Simpson's  rule. 


The  calculated  values  of  strain  by  Methods  A  and  B  for  the  recover},'  at 
decreasing  temperature  are  shown  in  Fig.  8.  Figure  8  shows  that  in  the  early 
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stage  of  recovery  the  test  data  lie  between  predictions  of  the  two  methods  but 
closest  to  Method  A.  However,  the  shape  is  predicted  more  accurately  by  Method  B. 

Conclus ions 

The  following  conclusions  may  be  drawn  from  the  experiments  and  analyses 
presented  in  this  paper. 

(1)  Results  of  the  constant  stress  creep  tests  at  different  temperatures 
indicate  that  creep  strains  are  separable  into  time- independent  and  time-dependent 
parts.  The  time- independent  part  of  the  strain  was  found  to  be  independent  of 
temperature  and  linearly  proportional  to  the  stress.  The  time-dependent  part  of 
the  strain  was  represented  by  a  power  law  of  time  with  a  constant  exponent  inde¬ 
pendent  of  temperature.  It  was  found  that  the  time-dependent  strain  was  a 
strongly  nonlinear  function  of  stress  and  temperature. 

(2)  Two  methods  were  considered  to  account  for  varying  temperature  by 
using  the  constitutive  relation  derived  from  different  constant  temperatures. 

One  method  employed  a  reduced  time  involving  a  shift  factor  which  was  a  function 
of  both  stress  and  temperature  history.  The  other  considered  the  strain  to  be  a 
function  only  of  the  current  temperature  not  its  past  history.  Both  employ  the 
modified  superposition  principle.  The  latter  method  agrees  somewhat  better  than 
the  former  with  the  trend  of  experiments  on  varying  temperature  and  recovery.  It 
is  also  the  easier  to  employ. 
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TABLE  I 

Tension  Creep  Experiments 


TEST 

TENSILE 
STRESS , 
ksi 

TEMPERATURE , 

°F 

e  =  e°  +  e+tn,  n  =  0.143 

o  . 

e  ,1 

+  g. 
t 

38 

1.0 

75 

0.195 

0.017 

40 

1.0 

100 

IT 

0.028 

51 

1.0 

130 

II 

0.053 

52 

1.0 

140 

II 

0.062 

53 

1.0 

150 

II 

0.068 

54 

1.0 

160 

M 

0.092 

63 

2.0 

75 

0.390 

0.052 

65 

2.0 

102 

II 

0.081 

41 

2.0 

115 

II 

0.100 

46 

2.0 

130 

II 

0.118 

72 

2.0 

135 

II 

0.125 

48 

2.0 

140 

II 

0.142 

77 

3.0 

75 

0.585 

0.083 

47 

3.0 

130 

II 

0.185 

50 

3.0 

140 

II 

0.267 

84 

4.0 

75 

0.780 

0.133 

92 

4.0 

102 

II 

0.178 

91 

4.0 

115 

It 

0.216 

78 

4.0 

126 

II 

0.268 

81 

4.0 

135 

It 

0.298 

97 

4.0 

145 

II 

0.328 

TABLE  II 

Repeated  Experiments  in  Tension 


TEST 

STRESS , 
psi 

TEMPERATURE  , 

°F 

STRAIN  AT  1  hr.. 
Per  Cent 

39 

2000 

75 

0.4341 

63 

ft 

II 

0 .4419* 

76 

If 

II 

0.4368 

41 

II 

115 

0.4902+ 

67 

tt 

II 

0.4925 

49 

II 

150 

0.6368 

114 

II 

II 

0.5894 

2 

3000 

75 

0.6439 

77 

II 

II 

0.6672+ 

87 

4000 

115 

1.0843 

91 

II 

II 

0.9957+ 

+  Used  in  Fig.  1,  2 
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Fig.  1.  Cree] 
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Fig.  3.  Creep  at  a  constant  stress  of  4  ksi  and  increasing  temperature 
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Fig.  6.  Vertical  shift  of  £+/tf  with  temperature  (see  Eq.  7d) 
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Fig.  7.  Recovery  at  zero  stress  and  at  the  indicated  temperature  following 
creep  for  one  hour  at  the  same  temperature  and  at  the  stress  shown 
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Fig.  8.  Recovery  at  zero  stress  during  decresing  temperature  following  creep 
at  a  constant  stress  of  4  ksi  during  increasing  temperature  (see  Fig.  3) 
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